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Single-sensor digital cameras capture image with the aid of masking the 
sensor surface along a colour filter array(CFA) such that every sensor pixel 
solely samples certain of three primary colour values i.e., R (red), 
G (green) and B (blue). To get a full-colour image, an interpolation method 
commonly referred in conformity with CFA demosaicking is required to 
estimate the other two contributions for producing a full-colour image. But, 
the clutter in imaging sensors not only corrupts the colour filter array but also 
introduces artifacts at some stage in the colour interpolation step and affects 
the characteristics of image. To acquire high quality full-colour image, a kind 
of viable and effective interpolation algorithm based over gradient is used. 
This technique can remove the noise effectively by retaining image border 
and detail data clearly. 
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1. INTRODUCTION 

With the evolution of science and technology in the defence and civil sectors, the colour camera 
with single CCD are extensively used as image input device. The colour image from the single CCD digital 
camera is referred to as the CFA (colour filter array) image. Currently, the colour image recovery algorithm 
primarily based on the CFA is widely used. The present writing put forward a lot over colour interpolation 
algorithms, namely adaptive interpolation technique [1], weight coefficient technique [2], interactive 
interpolation method [3], based on vector [4], and so on in an optimized way. However often used 
interpolation technique is bilinear interpolation technique [5] that belongs to the single channel independent 
interpolation method. In this technique the unknown colour factor among a point is computed generally by 
means of the average of adjacent same colour components. This approach runs faster, however ignores the 
detail data and the correlation between the three-colour channels, therefore the bilinear method frequently 
cannot achieve effective interpolation. Colour proportion constant method [6], it has an intense relation 
within different colour channels considering the correlation and the quality of the reconstructed image, was 
improved, but in fact still belongs to the class concerning bilinear method. The method primarily based on 
gradient [7], [8], researchers introduced the interpolation algorithm based totally on gradient, that can select 
the appropriate interpolation direction and can avoid the appearance of the zigzag pattern in the edge of 
recovered image. But, this approach does no longer consider the influence over noises of the image and 
accomplish the colour recovery distortion close by the noise. The proposed interpolation method can remove 
the clutter primarily based on gradient and effectively excerpt the impact of the noise by retaining the edge 
and the detail information of the image. 

Removal of noise in color image in an optimized way is achieved by using red component for the 
interpolation. Red component is used for interpolation from the R, G, and B components for removing the 
noise in order to optimize the computations required. As R and B components are accounted each as 4 of the 
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total number of the pixels in Bayer CFA pattern. Whereas the G component is 1⁄2 of the total number of 
pixels. Using G component for the interpolation, in order to remove noise takes more time as it includes more 
computations compared to R and B component. 


2. CFA IMAGE COLOR RECOVERY METHOD 
2.1. CFA (color filter array) image 

There is only colour component gray value on each lattice point in the CFA (colour filter array) 
image. Because of the human eye photosensitive characteristic, at present the GRGB colour swatches is most 
commonly used, namely Bayer colour filter array, as shown in Figure 1. 


Figure 1. Bayer CFA pattern 


It uses a group of red and green filter or a group of blue and green filter by turns to obtain image, the 
number of green pixels are half part over the other pixels, and the red and blue then each for 1/4. Due to the 
green component accounted for half of the total, hence it has more detail information over image, therefore, 
the interpolation algorithm begins mostly advance from restoring G component. 


2.2. Filter out noise method based on gradient 

The technique based on gradient, does not consider the impact of noise to algorithm, then the image 
entails G11 noise, if the clutter as colour information involved in calculation after recovering image, not only 
makes the colour distortion, but also using the information of four restore point close to the noise, their 
colour component also can appear distortion. Hypothesis, Gi,j, is a high frequency clutter point, then, G 1-i,j 
Gl+i,j, Gij-1 Gij+1 and Gi, , their G factor will appear distortion. Therefore, it is important to remove 
clutter for getting better colour image, but generally the median and mean filter is used for the gray image 
method, are not appropriate for CFA distribute images. This paper study two techniques primarily based on 
gradient form Hibbard [3] and Laroche [4], through the gradient of the calculation results in both technique 
connected, eliminating the impact of the noise. 

As shown in Figure 1, B;j, is a point of B component in the image, in order to restore Gj , Gj, says 
the value of G component in this point. In (1) Al is horizontal internal gradient and B1 is vertical internal 
gradient, through calculate one order differential such as formula (1). 


Ay = |Gij-1 — Gi j+ | (1) 
B, = [Gi-ty = Gi+1,j | 


In (2) A2 is horizontal external gradient, B2 is vertical external gradient, through calculation two order 
differential such as formula (2). 


l Az = |2 X Bij — Bi j-2 — Bi j+2l 2) 
B, = |2 X Bij — Bi-2,j — Bisa, 


According to the gradient results of internal and external two layers, to locate edge information of 
image if really exist, or have the influence of the noise point. Set TH is enumeration variable, for being 
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clutter that in the up and down or so said at the point of G,, and for having horizontal or vertical edge 
information in this point. To compute the formula concerning TH. 


up (Ay < By) & (Az = By) & (max|2 x Giy — Gi-1j — Gira j| =j — 1) 
down (A, < B1) & (Az = B2) & (max|2 x Giy — Gi-sy — Gim j| =j +1) 
left (Ay > By) & (Ap = By) & (max|2 x Gy j — Gij-1 — Gij41| =i-1) 
TH =} right (A, > By) & (Az = B2) & (max|2 x Gy; — Gi j-1 — Gi jai] = i+ 1) (3) 
no (Ay = By) & (Ap = B2) 
level (A, > B1) & (Az > Bz) 
erect (A, < B1) & (Az < B3) 
((=2,4,6...57=13,5...3x=i-1Lit+ly=j-1j+1 


In (3) max|2 x Gj, y-Gi-1, ;-Gisi, j |=j+1 is the column location of point, whose distance is farthest 
between G;.;,; and G;,,;, max|2 x Gx, j-G;, ;.1-Gi, j+1 | whose distance is farthest between, Gij-ı and Gij1, If TH 
equals up then it indicates , Gjj.; is the noise, down indicates , Gij+ is the noise, left indicates G;.,;, is the 
noise right indicates Gj,;; is the noise, no indicates that there is no noise and no edge, level indicates that 
there is an edge in vertical direction, erect indicates that there is an edge in horizontal direction. 

So the finally calculate formula of G;;, like formula (4) below. 


(Gi-a,j + Gizaj + A TH = up 
(Gi-1 + Gisaj + Gij-)/, TH = down 
(Gisa j + Gij-1 + A TH = left 
p (Gia + Gij-1 + A TH = right (4) 


(Gizi, + Ginny + Gij-1 + Gijs)/ TH = no 
A = 


(Gia; + Sns), 


TH = level 


(Gij-1 + Gisj+1)/ 
2 TH = erect 


The noise removal after the G;j restored is as follows, 


Gij-15 Gij TH=up 

Gi ja = Gij TH = down 5 
Gi-1,j = Gij TH = left 

Gi+aj = Gij TH = right 


3. RESULTS 

The test image which is used to apply the denoising technique is of the size 925(H) x 590(V). The 
Figure 2. is the original CFA image with noise and Figure 3. is the denoised image effectively filter the noise 
and make the image look more refined. 

The quality of two images is measured using SNR. Here Figure 2 is the image before applying the 
filter technique and Figure 3 is the image after applying the filtering technique. 
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Figure 2. CFA image with noise Figure 3. Denoised image after applying the 
algorithm 


SNR for a given image can be computed using the expression, 


SNR=/VLSD? (6) 


Where u is average gray of colour image and LSD max is local variance maximum. 


4. 


Table 1 shows the SNR of two Figures. 


Table 1. SNR of two Figures 


Figure SNR (existing) SNR (proposed) 
2 17.09333 33.165210 
3 21.66333 39.215886 
CONCLUSION 


The CFA image colour interpolation method introduced in this paper used filter out noise 


interpolation method based on gradient in an optimized way to avoid noise on the colour recovery influence. 
This method has wide application in defence and civil sectors which improves the signal to noise ratio of a 
colour image and has a wide application prospect. 
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